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Abstract 

In this article an upscaled model is presented, for complex networks with highly clustered regions exchanging 
some abstract quantities in both, microscale and macroscale level. Such an intricate system is approximated 
by a partitioned open map in IR 2 or IR 2 . The behavior of the quantities is modeled as flowing in the map 
constructed and thus it is subject to be described by partial differential equations. We follow this approach 
using the Darcy Porous Media, saturated fluid flow model in mixed variational formulation. 
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1. Introduction 

A highly clustered network, is a graph such that its clustering coefficient is close to one. The clustering 
coefficient [T] of a vertex v in a graph is defined as the number of triangles connected to v divided by 
the number of triples where v is incident on two edges (triples “centered” at v). The network clustering 
coefficient is the average value of the clustering for every node in the graph. In this work we consider a 
simple, connected network where two nodes are linked if and only if they exchange some abstract quantities 
and containing highly clustered regions or universes (communities, depending on the context). Hence, these 
universes also exchange the quantities of interest with other clustered regions in the network or multiverse, 
but at a different scale. Therefore, a macroscale graph representing each clustered region with a node and 
connecting two nodes whenever the universes they represent trade these quantities, is a natural upscaled 
model. Additionally, the number of nodes in each universe is very large compared to the number of clustered 
regions in the multiverse and the exchange of abstract quantities between universes is very large compared 
to the exchange amongst the individual nodes in the universe. A natural example is the exchange of goods 
and services amongst members of a nation and between the nations; other real world systems resembling 
these characteristics can be found in biotechnology [2], social and economic networks [3], the Internet Pi, 
etc. The aim of this article is to provide an overall description of how these quantities are exchanged at 
the macroscale level. Such necessity has already been stated implicitly in [3] according to the quote: “... 
what makes these networks complex is that they are generally so huge that it is impossible to understand or 
predict their overall behavior by looking into the behavior of individual nodes or links...”. On the other hand, 
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a highly descriptive model becomes impractical for complex networks because of the elevated computational 
costs, numerical instability and low quality solutions introduced by large-scale computations. 

The relationship between PDE and graphs has been subject of study in recent years. Most of the work 
has been done to provide the basic definitions of the domain associated with the graph and the strong 
differential operators defined by the PDE of interest, see 0 for a global survey on the field. Authors 
commonly choose the 1-D simplicial complex given by an embedding in IR, V of the studied graph, define 
strong operators in the edges and matching conditions on the vertices, together with appropriate function 
spaces. The mathematical approach is mainly classical and the technique heavily relies on the eigenfunction- 
eigenvalue expansion methods and/or maximum-minimum principles; as an example of this treatment see 
|Hj, for a broad exposition see [7]. The results depend crucially on the geometry of the embedding however, 
it is not clear how to make such choice. Seeking to gain independence from this limitation another approach 
consists in defining discrete difference operators, mimicking the properties of the PDE operators. Again, 
the subsequent mathematical treatment depends on eigenvector-eigenvalue methods and their properties, 
see [§] for a deep discussion, followed by the construction of Green’s functions, see [9] as an example. Yet 
an intermediate approach addresses time-evolution problems using discrete models for space, as imposed by 
the graph itself, and continuous evolution in time under the hypothesis that the underlying combinatorial 
structure of the graph remains stationary, see 0 for this view. In contrast with the previous achievements 
the present work preserves the continuous definitions for the operators in the PDE, and adapt the domain 
associated with the graph allowing for the use of weak variational formulation results. In that sense this 
article provides a dual approach to the previous results, however our main motivation is to attain upscaling 
criteria for highly complex networks which is a strong necessity as previously discussed. Next, we describe 
the model introduced in this paper. 

Each clustered region will be modeled by an open bounded, simply connected set in IR. v where every 
point will represent an individual/molecule. In order to approximate the clustering, we propose that every 
point will exchange these quantities with every element in a small neighborhood. Following this concept 
if two universes exchange the quantities in question then, the sets representing them will share a non- 
negligible boundary. Moreover, according to [Jj the highly clustered regions are the ideal medium for rapid 
communication between the nodes. Therefore, we propose that the observed quantities can be realized as a 
fluid flow phenomenon in the modeling open set. For simplicity we adopt the stationary, saturated, Darcy 
Flow model 0 to approximate the exchange of the aforementioned quantities within the region 

au+Vp + g = 0, (la) 


V-u-F = 0 infi. (lb) 

Here, a is a positive coefficient describing the resistance to the flow of the medium. Additionally, coupling 
conditions will be introduced to describe the exchange between universes, as well as boundary conditions for 
the multiverse overall behavior. This is a PDE problem defined on a domain associated to a graph, however 
in order to model successfully the exchange of two quantities simultaneously, the available tools of analysis 
m demand the underlying graph to be bipartite; this additional hypothesis will be necessary and included 
in Section |4j Finally, the coefficient a in (la) can be interpreted as resistance to the flow within a network 
e.g. fees and taxes slowing down the exchange of goods and services in a nation state, paradigms impeding 
to permeate new ideas in a social network [3] , band width limiting the diffusion of information through the 
Internet [5], etc. 

The paper is organized as follows. In Section[2]we list the results and concepts needed for the exposition. 
Section [3] defines the types of domain to be associated to the graph and proves their existence. Section 
[4] introduces the PDE model together with the necessary geometric associated notions, it also shows the 
formulation of the problem, proves its well-posedness and recovers the strong form. Finally, Section [5] 
presents the final discussion, and future work. 
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2. Preliminaries from Graph Theory and PDE 

2.1. Preliminaries from Graph Theory 

We begin this section with the basic, necessary definitions from graph theory El- 
Definition 1. Let G = {V. E) be a graph 

(i) The degree of a vertex v, is the number of edges that have an endpoint at v. 

(ii) A walk in G from vertex vo to vertex v 3 is an alternating sequence 

(co, ei, v\ , e*2 , • * ■, v n —i , e n , u n ), 

of vertices and edges such that the endpoints of the edge e, are Vi-\ and V{ for all i = 1,2,..., n. 

(iii) A path is a walk with no repeated edges and no repeated vertices, except possibly the initial and final 
vertices. 

(iv) A walk or path is trivial if it has only one vertex and no edges. 

(v) A cycle is a non-trivial closed path i.e. it starts and ends on the same vertex. 

Definition 2. A self-loop is an edge that joins a single vertex with itself. A multi-edge is a collection of 
two or more edges joining identical vertices. A simple graph has neither self-loops nor multi-edges. 

Definition 3. (i) A graph is connected if for every pair of vertices u and v there is a walk from u to v. 

(ii) An edge e is a bridge, if the graph G — e is not connected. 

Definition 4. A graph with no cycles is a forest, if additionally the graph is connected it is said to be a 
tree. In a tree, a vertex of degree one is said to be a leaf. 

The following is a well-known result about trees, E- 
Proposition 1. A tree with at least one edge has at least two leaves. 

Definition 5. A bipartite graph is a graph whose vertex set V can be partitioned in two subsets U, W 
such that each edge of G has one endpoint in U and one endpoint in W. The pair U, W is called a (vertex) 
bipartition of G', and U and W are called the bipartition subsets. 

Next we recall several definitions about planar graphs E- 

Definition 6. (i) A graph is said to be embeddable in the plane, or planar, if it can be drawn in the 

plane so that its edges intersect only at their ends. 

(ii) A planar embedding of a graph will be referred to as a plane graph. 

(iii) A plane graph G partitions the rest of the plane into a number of arcwise-connected open sets. These 
sets are said to be the faces fo G. 

(iv) We say that a vertex v of a plane graph G is an outer vertex, if it belongs to the boundary of the 
outer face of G. 

Definition 7. Let G be a plane graph, 

(i) The dual graph G* is defined as follows. Corresponding to each face / of G there is a vertex f* of 
G* and corresponding to each edge e of G there is an edge e* of G*. Two vertices f* and g* are joined 
by the edge e* in G* if and only if their corresponding faces / and g are separated by the edge e in G. 
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(ii) The plane dual of the plane graph G is a natural embedding of G* in the plane. It is obtained by 
placing a vertex f* in the corresponding face / of G, and then drawing an edge e* in such a way that 
it crosses the corresponding edge e of G exactly once and crosses no other edge of G. We refer to such 
a drawing as plane dual of the plane graph G. 

Definition 8. A curve in ]R 2 or IR 3 is the continuous image of a closed interval, we say that a curve is 
simple if it does not intersect itself. 

We close this section recalling two well-known results GH 113- 

Theorem 2. A plane graph G is connected if and only if it is isomorphic to its double dual G**. 
Theorem 3. Every finite graph is embeddable in 1R 3 . 

2.2. Preliminaries from PDE 

We start this section introducing the general notation. In the present work vectors are denoted by 
boldface letters as are vector-valued functions and corresponding function spaces. The symbols V and V- 
represent the gradient and divergence operators respectively. The dimension is indicated by N which will 
be equal to 2 or 3 depending on the context. Given a function / : IR. Y —> IR, then f M f dS denotes the 
integral on the N — 1 dimensional manifold M. C ]R, Y . Analogously, j 4 /dx stands for the integral in the set 
A C IR/^; whenever the context is clear we simply write f A f. The symbol V denotes the outwards normal 
vector on the boundary of a given domain O C IR, Y . Given an open set M of IR/^, the symbols || • ||o,m ; 
II • IIi.m, || • Hi/ 2 ,dM, II • 11-1/2, 3 M and || ■ || Hdlv (M) denote the L 2 (M), H l {M), H 1 / 2 {dM ), and 

Hdiv(Af) norms respectively, while \M\ represents the Lebesgue measure of M in IR 2 or IR, 3 depending on 
the context. 

Next, we present the general abstract problem to be studied in this article. Let X and Y be Hilbert 
spaces and let A : X —> X', B : X —> Y' and C : Y —5- Y' be continuous linear operators, we are to study 
the following problem 


Find a pair (x, y) € X x Y : Ax + B'y = F\ in X', 

-Bx + Cy = F 2 inY'. 


( 2 ) 


Here F\ £ X' and F 2 € Y'. Several variations of systems such as the above have been extensively studied, 
we present below a well-known result m to be used in this work. 


Theorem 4. Assume that the linear operators A : X —> X', B : X —> Y', C : Y Y' are continuous and 


(i) A is non-negative and X.-coercive on ker (B). 


(ii) B satisfies the inf-sup condition 


inf sup ■ 
ye Y x gx 


|g*(y)| 

l x l!x llyllv 


> 0 . 


(3) 


(in) C is non-negative symmetric. 


Then for every F\ £ X' and F% £ Y' the problem ([2]) has a unique solution in (x, y) £ X x Y, which satisfies 
the estimate 


ll x llx + l|yII y < c(||Fi||x' + II^Hy')- 


(4) 
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(a) Downscaling Map. 


(b) Tubular e-Map. 


Figure 1: Figure (a) depicts a downscaling map example for a given plane graph represented in dotted line. Figure (b) depicts 
a tubular map example for the same given plane graph. 


3. The Graph Domain 

This section is aimed to the construction of a particular topological domain for a given plane graph. The 
domain must be suitable for setting a PDE problem. To this end, we introduce two paramount definitions 
of domains associated to graphs depicted below 

Definition 9. Let G = (V. E) be a connected graph embedded in either ]R 2 or ]R 3 such that its edges 
are simple curves. Denote by v, e the points and lines representing the vertices and edges of G and 
G '= f (U{£ : v € V},\J{e : e G E}). Let e > 0 be such that the collection of balls {B(v,e) : v & V} 
is pairwise disjoint. Define 

(i) The tubular e-region 

Uq d = f {x : d(x, G) < e}. (5) 

(ii) For each edge e G E choose a smooth simple curve e', approximating e. Denote by l e a secant line (or 
secant plane) through the midpoint of e' and let C e be the connected component of Uq fl £ e containing 
such midpoint. 

(iii) Let {C e : e G E} be as defined above. For each v G V and w adjacent to v, let U r u w be the tubular 

e-region corresponding to the induced subgraph G VtW *= ({v,w},vw). The set C vw divides U{, w into 
two open regions, one containing v and one containing w denoted by H(v, w) and H(w, v) respectively. 
Define the starred region of v by 

W V A = (J H(v,w). (6) 

w 6 V: vw 6 E 

(iv) The collection { U,l : v G V} is said to be a tubular e-map, or simply a tubular map, of the graph G. 
Now, we introduce the concept of downscaling map. 

Definition 10. Let G = (V. E) be a plane connected graph, wc will say that a downscaling map of G is 
a collection of bounded open sets {O v : v G V} called regions, such that 

(i) v G O v for all v G V. 
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(ii) If v w then O v fl O w = 0. 

(iii) O v is simply connected for all v £ V. 

(iv) The set O defined by 

O d =int (cl IJ oA (7) 

V vev ) 

is simply connected. We define O as the domain of the downscaling map. 

(v) Two elements of the collection share non-negligible boundary if and only if the vertices they contain 
are connected in the graph, i.e. \dO v H dO w \ > 0 if and only if vw £ E. 

(vi) If v is an outer vertex then | dO v fl dO\ > 0. 

Finally, we will say that the regularity of the map is given by the lowest degree of regularity of its elements. 
Remark 1. Notice the following 

(i) A tubular map satisfies all the conditions of a downscaling map, except possibly for the global simply 
connectedness condition (Definition |To]Jiv| ) ). 

(ii) A tubular map of a plane graph defines a downscaling map if and only if the graph is a tree. 

The next two results are central in proving the existence of a downscaling map for a simple, plane, connected 
graph. The intuitive idea and technique are depicted in Figure 1(b)- 

Lemma 5. Let G = (V) E) be a connected, simple, plane graph such that no bridges are in the boundary of 
its outer face. Then, there exists a downscaling map for G. Moreover, this existence can be attained for any 
chosen level of regularity. 

PROOF. If G has no edges then it must be a single vertex v, thus an open ball centered at v will satisfy the 
definition of downscaling map. We will henceforth assume that G has at least one edge. 

Let G* be the plane dual graph of G drawn with disjoint simple curves as edges. The regions defined by 
the faces f* of G* are the natural candidate to define a downscaling map of G , however they fail because of 
two reasons. On one hand, according to Theorem [2] we know that the double dual G** is isomorphic to G. 
In particular, the outer face of G* must contain a vertex of the plane graph G. On the other hand, given a 
face f* of the plane dual containing an outer vertex v, it would not necessarily hold that \df* 0 <9/q | > 0, 
where /g is the outer face of G*, as demands the condition ( jvij ) in Definition 

We overcome the first deficiency as follows, let vq be the unique vertex in G contained in the outer face 
/g of G* and let C be the cycle in G bounding its outer face /o. Clearly vq belongs to C , now let e > 0 
be such that the tubular e-region of the induced graph C — vq is completely contained in (/g) c . Define the 
region 

C’o d = f fo n U B ^ e ^ 

x e(/o) c 

Since /o is the outer face, it is clear that /q and (J x £ B(x , e) is simply connected. Therefore Oq is open 
bounded, simply connected and contains Vo- For the second deficiency, let xq be the point representing the 
outer face fo in the plane dual G*. Now let S > 0 be such that cl B(xq,S) C fo, then the collection 

{/* - cl B(x 0 ,S) : f* face of G*} U {O 0 — c\B{x 0 ,5)}, 

constitutes a downscaling map for the graph G. 

Finally, since the smooth curves are dense in the plane, it is clear that the boundaries of the elements 
of this downscaling map can be continuously deformed, to define a new downscaling map with any required 
level of regularity. □ 
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G H, 



(a) Starting Graph G. 



(b) Dual Graph of H = G(V, E - F). 


Figure 2: Figure (a) depicts a starting graph G and its set of outer bridges F, it also shows in dashed line the tree Hi removed 
from G, as in the proof of Theorem [7] Figure (b) depicts the construction seen in Lemma [ 5 ] for the graph H depicted in dotted 
line. The construction of the domain Oo and the removal of the ball B(x 0 ,5) are depicted in dashed dotted line. 


Proposition 6. Let G be a connected, simple, plane graph, let F be the set of bridges in the outer face of 
G, and Hi,...,Hk be the connected components of the graph defined by removing the edges in F. Let G be 
the graph whose vertex set is {rq, ...,Vk}, where Vi is connected to Vj if and only if Hi is connected to Hj by 
an edge in F. Then, G is a tree. 

PROOF. If G contains a cycle, the removal of one edge in such cycle does not disconnect G. Hence such 
edge is not a bridge contradicting the definition of F. □ 

Finally, we present the main result of this section. 

Theorem 7. Let G = ( V , E ) be a connected, simple, plane graph. Then, there exists a downscaling map 
with smooth boundaries for G. 

PROOF. First, if G has no cycles then it is a tree and its tubular map constitutes a downscaling map. Hence, 
from now on we assume that G has at least one cycle. 

Let Hi,...,Hk be the components and G be its associated graph as defined in Proposition [6] We will 
proceed by induction over k. If k = 1, then Lemma [5] provides the required downscaling map. If k > 1, 
then renumbering the components if necessary, suppose that Hk is such that Vk in G is a leaf. Let A be 
a downscaling map with smooth boundaries for Hk and let M be its domain. Analogously, let 0 be a 
downscaling map with smooth boundaries for the graph such that includes H 1 ,..., Hk- 1 together with every 
bridge connecting these components and let O be its domain. By deforming their boundaries, it can be 
additionally assumed that d(0,M) > 0, while preserving the smoothness of their boundaries. 

Let v and u be the vertices such that vu is the unique bridge connecting Hk with some component Hi 
with 1 < i < k. Clearly v € Hk and u € H^. Denote by O u and M v the regions of 0 and A containing u and 
v respectively. Consider the graph G consisting of the vertices u, v and the edge uv. Take two points p and 
q in the interior of the one dimensional manifolds 80 H dO u and dM fl dM „ respectively. Choose a simple 
curve C connecting u and v, which intersects dO and dM at the unique points p, q and such such that 
d[C , O') and d(C,M') are positive; where O' d = f O — O v and M' *= f M — M v . Choose e > 0 strictly less 
than min{d(C, 0'),d(C, M')} and such that if Q denotes the e-map of G, then it verifies that Q fl O C O u 
and Q n M C M u - Redefine O u , M v in order to include the regions in the tubular e-map Q, corresponding 
to u and v respectively. If necessary, deform the boundaries of O u and M v to attain the required smoothness 
and denote the outcome by O u and M v respectively. The collection A — {AIl,}, 0 — {O u } together with 
{O u } and {M v } constitutes a downscaling map for G with the required regularity. □ 
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4. The Downscaled Bipartite Model 


4-1. Geometric Setting and Modeling Function Spaces 

In this section we give the geometric setting for the formulation of the problem. 


Definition 11. Let 12 be a connected bounded region with smooth boundary, let Q = {K : K £ Q} be a 
bipartite map and denote by Q\ = {L : L £ Q i}, Qn = {M : M £ C/ 2 } the bipartition, or bi-coloring, of the 
map. 

(i) For each region K £ Q denote by v the outer normal vector to its boundary dK. 

(ii) For each region K £ Q denote by n : dK —> IR lV the “normal” vector by 

(v(x) I\ £ Gi, 

n(x) d = f < —V(x) I\ £ Q 2 and x £ dK n 12, (8) 

[ v(x) I\ £ O 2 and x £ dK n 512 . 

(hi) Define 12i = f U(L : L £ £ 1 }, D 2 d = (J {M : M £ Q 2 } . 

(iv) Define T d = f |J{5A' : K £ Q} - dfl. 

(v) Define 


£ = f \J{dL fl dM : L £ Qi, M £ Q 2 } 


{0 : a is the interface of two regions of different type} 
{cr : cr defines an edge in the graph G of the map Q}. 


(9) 


Remark 2. Notice that 


(i) Simple connectedness upon the domain is not required in order to include the tubular map. 

(ii) In agreement with the previous section notice that a bipartite map will always be induced by a bipartite 
simple graph. 


In order to successfully associate a well-posed problem we endow the Problem 0 with boundary conditions 
(10c), (lOf), together with the exchange interface conditions of normal flux balance (lOh) and normal stress 
balance (lOh). This gives the following strong problem 


«u-| + Vpi + g = 0, 

(10a) 

V • u 2 = F in 12 2 . 

(lOe) 

V • Ui = F in 12 1 . 

(10b) 

u 2 • n = 0 on 512 2 D 512, 

(lOf) 

Pi = 0 on 512i n 512. 

(10c) 

ui • n - u 2 • n = /3 p 2 + fa , 

(10g) 

a u 2 + Vp 2 + g = 0 , 

(lOd) 

P2 ~ Pi = fs on T . 

(lOh) 


Hypothesis 1. It will be assumed that the storage exchange and the friction coefficients, j3 : F —>• [0, 00 ), 
a : 12 —> (0, 00 ) respectively, verify that /? £ L°°(T ), ||/3 l r || L i( r ) > 0 and a £ L°°( 12), ||^||z,°°(n) > 0. 








In order to introduce the modeling spaces used in the formulation we first notice that {L : L £ Q 1 } and 
{M : M £ Q 2 } are the simply connected components of Sli and fl 2 respectively. Then, 

Hm^) = 0 H div (L), H 1 ^) = 0 H\M). 

l e Qi Me 62 

The following space is introduced in order to couple adequately, the action of the pressures traces in the 
variational formulation 

E(Cl 2 ) d = f {<? G H 1 ^) : qtdMndL G H 1/2 (dL) for all (L, M) eftx Q 2 ) 

= {q£H\n 2 ) :ql r £ H 1 / 2 ^)}. 

We endow E(Ll 2 ) with the H 1 ^ 2 ) inner product. It is direct to see that E{Ll 2 ) is a closed subspace of 
H 1 ^ 2 ) and consequently a Hilbert space. Now define 

V(fi 2 ) d = {v £ L 2 (H 2 ) : v 2 = Vq 2 for some q 2 £ E{Q 2 )} = V(£(H 2 )), (12) 

endowed with the L 2 (fl2) inner product. Next we show a necessary result. 

Lemma 8. Let E(Tl 2 ) and V(fl 2 ) be as defined in respectively, and define 

Efifih) d = {q 2 £ E(n 2 ) : jf q 2 = o}. (13) 

Then, 

(i) There exists a constant C > 0 depending only on the domain Ll 2 such that 

IM|i,n 2 < C'||Vr 2 ||o,n 2) for all r 2 £ H. (14) 


(ii) The space V(fl 2 ) is Hilbert. 


PROOF. (i) Clearly E 0 (n 2 ) is closed and because of the matching property of traces for elements of E[fil 2 ) 
the application r 2 1 —»• || V ?~2 ||o,o 2 is norm in Eq{L1 2 ). Due to the Rellich-Kondrachov Theorem this norm 
is equivalent to the standard one in Eq(L1 2 ) i.e., there exists C > 0 depending only on the domain fl 2 
satisfying the statement (14). 


(ii) Clearly, it is only necessary to check that V(f2 2 ) is complete. Let {vj : n £ IN} be a Cauchy 
sequence in V(D 2 ), then there exists a sequence {pVq : n £ IN} in E(Ll 2 ) such that Vq% = vlf, then, 

the function rV; qff — f n ^ q 2 belongs to £ 0 (^ 2 ) and Vr 2 = v 2 . Due to the previous part it 
follows that the sequence {r 2 : n £ IN’} C Eq(L1 2 ) is Cauchy, consequently it converges to an element 
r 2 £ E 0 (n 2 ) C £(fi 2 ). Finally, since the gradient map V from £(fl 2 ) onto V(fl 2 ) is continuous, the 
result follows. □ 


Now we are ready to introduce the functional setting of the problem, we have 


X d =H div (0 1 ) xE(Si 2 ), 

(15a) 

Y d = V(D 2 ) x L 2 (fl 1 ). 

Endowed with their natural inner product and norms 

(15b) 

[vi,g 2 ] x = f {l|vi||^ dlv(0l) + I 92 ffi ( n 2 )} 2 : 

(15c) 

[ v 2,9i] Y = f { V2 |L2 ( n 2 ) + Mli* ( n x )} 2 - 

(15d) 
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Remark 3. (i) Notice that the definition of spaces gathers the functions of high regularity in X and the 

functions of low regularity in Y. This choice is made on one hand to satisfy the hypotheses of Theorem 
[4]and, on the other hand, to preserve the remarkable aspect that the underlying modeling spaces X and 
Y are free of coupling conditions. This approach will lead to a version of mixed formulation different 
from the one presented in m and nn. which shares the coupling-free spaces feature. 

In order to avoid heavy notation, in the sequel we adopt the following conventions 


(ii) Let A be an open bounded set, v £ Hdi v (A) and q £ H 1 ( A), then we denote 


[ ( v-n)qdS d = (v • n, q) 

Jd A 


/ 2 (0A),ff-l/2(dA)' 


(iii) Since T = (J ff _ £ a we denote 


J (vi • n) q 2 dS d = f ^ J (vi • n) q 2 dS. 


(16) 


(17) 


<7 £ £ ' 

4-2. Weak Formulation of and Well-Posedness of the Problem 


In this section we present a particular mixed-mixed formulation for the problem (10). 


Find ([ui,p 2 ], [u 2 ,Pi]) £ X x Y : / au r vi+ / f3p 2 q 2 dS - / (ui • n) q- 2 dS + / p 2 (v 1 -n)dS 

Jr 2 X Jr Jr Jr 

- [ Pi v • Vi - [ U 2 • V <?2 = f F q 2 — f g ■ Vi + [ /s (vi • n) dS - [ fa q 2 dS, (18a) 
Js 2 i Jn 2 Jn 2 Jn i Jr Jr 


/ V • Ui (7i + / Vp 2 -v 2 -|-/ a u 2 • v 2 = / Fq!- g-v 2 for all ([vi, q 2 \, [v 2 , gi]) £ X x Y. (18b) 
JJ f22 J" ^1 ” ^2 

Define the operators A : X —> X', B : X —> Y' and C : Y —► Y' by 

-4[vi,g 2 ], ([wi,r 2 ]) d = f [ a vi • wy + f fiq 2 r 2 dS- I (vi -fi)r 2 dS+ f q 2 (wi • n) dS, (19a) 
J Oi Jr Jr Jr 

B[vi, q 2 \, ([w 2 , ri]) d = f f V • v x r x + f V q 2 ■ w 2 , (19b) 

J J ^2 

C[v 2 , gi]([w 2 , n]) d = f [ av 2 ■ w 2 . (19c) 

Jih 


Thus, the Problem (181 is equivalent to 


Find a pair ([u 2 , pi], [u 1( p 2 ]) € X x Y : A[u- 2 , px] + p 2 ] = F 1 in X', 


( 20 ) 


-B[u 2 , pi] +C[ui, p 2 ] = F 2 inY'. 

Where F\ £ X' and F 2 £ Y' are the functionals defined by the right hand side of ( |18a ) and (18b) respectively. 
4-2.1. Inf-Sup Condition of the Operator B and Coerciveness of the Operator A on X n ker(B) 


Lemma 9. The operator B : X —> Y' defined in equation (19b) is continuous and satisfies the inf — sup 


condition i.e., there exists a constant C > 0 depending only on the map Q such that for every [w 2 ,ri] £ Y 
there exists [vi, q 2 ] £ X satisfying 


B [vi, q 2 \([w 2 ,ri\) > C || [v 1; q 2 \ || x || [w 2 , ri]|| Y . 
Moreover, the constant C > 0 is independent from [w 2 ,ri]. 


( 21 ) 
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PROOF. It is direct to see that the operator B is continuous. Now fix [w 2 , r\] £ Y, for each polygon L £ Qi 
let £ Hq ( L ) be the unique solution of the local homogeneous Dirichlet problem 

- V • V£l = r\ 1 L in L. Cl = 0 on dL. (22) 


Taking vl = VCl due to Poincare inequality we observe that ||vi|| Hdiv (L) < Cl ||»t1l||o,l where Cl 
depends only on the diameter of the simply connected region L. Therefore, the function vj d = f J2l e Gi v l -"-l 
clearly belongs to H div (fii) and || Vj. || Hdiv(ni ) < (max Le6l C L ) ||r-i||o ifil . 

Let w 2 £ V(H 2 ), by dehnition there must exist r) £ E(PL 2 ) such that V 77 = v 2 . Then, q 2 d = f 77 — jjyn/p 2 V 
belongs to the space Eq(£1 2 ) (defined in (fl3])) , it satisfies that Vr 2 = w 2 and due to the Inequality ( fl4| ) , it 
holds that ||<? 2 1|i,r 2 2 < C ||v 2 || 0 ,o 2 with C > 0 depending only on the domain 0 2 . 

Hence, the pair [vi,g 2 ] belongs to X and C || [ v o 92 ] || x < ||[w 2 ,n]|| Y for an adequate constant C > 0 
independent from [w 2 ,ri]. Moreover, 


H[vi,g 2 ]([w 2 ,n]) = [Hw 2 ,n]|| > C ||[v 1 ,g 2 ]|| x ||[w 2 ,n] 


l Y ' 


The inequality above yields the inf-sup condition (21). 


□ 


Proposition 10. The operator A : X —► X' defined by (19a) is X.-coercive on X fl ker(S) i.e., 
A[vi 1 q 2 )([v 1 ,q 2 }) > C||[vi,g 2 ]||^., for all [vi,q 2 ] £ Xnker(H). 


(23) 


Where C > 0 is an adequate constant depending only on the map Q and the storage coefficient (3. 


PROOF. The continuity of the operator A follows applying the Cauchy-Schwartz inequality on each of its 
summands and noticing that the boundary terms involved can be controlled by the norm || • |j x . For the 
coerciveness of the operator, let [v-,, q 2 ] £ X D ker(H) then 


S[v 1 ,q , 2 ]([w 2 ,ri]) = 0 for all [w 2 ,ri] £ Y. 


(24) 


In particular, testing (241 with [0,ri] £ Y we conclude that V ■ Vi = 0 since ri is an arbitrary element in 
L 2 (H 1 ). On the other hand, clearly Vq 2 £ V(fl 2 ) and the pair [Vg 2 , 0] £ Y is eligible for testing (24). The 
test yields V<7 2 = 0 i.e., q 2 is constant inside fi 2 . Hence 


Jet = ll/il g(‘ (rl tefc, = ll«i « 112 


|h 2 


2 II l.Qo 


Using the previous observations we get that 

-4[vi,<72]([vi, 9 2 ]) = / a Vl -vi+ / /3q%dS 
Jo. 1 Jr 


> 


1 -1 

a 


„ ,,2 , 11/53Lr||Li(r) ||2 

L~ ( n) Vl Hdlv ( ni ) + \Sh\ ll92,ll ’ n2 


where C = min (11- 


> C'||[ V 1:92]|| x , 

|fl 2 1 — 1 1|/3 lr||Li(r) }■ This completes the proof. 


111—1 


□ 


Theorem 11. The Problem (|20j) is well-posed and there exist a constant C > 0 depending only on the map 
Q such that 

11 [ u 2) Pi]|| x + ||[ui, p 2 ]|| Y < C'(||Fi|| X ' + ||.F 2 ||y')- 


Proof. The proof is a direct application of Theorem [4] as all the required hypotheses are satisfied. □ 
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4-3. Recovering the Strong Problem 

We begin this section with the strong problem that is modeled by the weak variational formulation (18). 
The process shows mild restrictions on the forcing terms which are characterized below. 


Theorem 12. The solution of the weak variational problem (18) is a strong solution of the problem (10) with 


the forcing gravitation term g in the equation (10a I replaced by P g, which denotes its orthogonal projection 


onto the space V(n 2 ). In particular if gln 2 £ V(f2 2 ) the weak solution is exactly the strong solution. 


and 


PROOF. First we focus on recovering the constitutive and conservative equations (10a), (10b), (|10eJ) 
(lOd), on the domains fli and 112 respectively. 

Choose [0,<7i] £ Y and test equation (18b) to get f n V • Ui gi = f n F q-\. for all qi £ L 2 {Vl i). This 
yields the strong conservation equation (10b). Next choose [v 2 ,0] £ Y and test equation (18b) to get 
/ 02 (Vp 2 + au 2 ) • v 2 = — g • v 2 . Clearly Vp 2 + au 2 £ V(fl 2 ) and the previous equality holds for all 
V 2 £ V(n 2 )- From here it follows that Vp 2 + au 2 = — Pg, where P g is the orthogonal projection of g onto 
V(fl 2 )- This gives the constitutive Darcy equation (lOd) with the forcing term g replaced by its projection 
Pg. Now let <f> £ then, testing equation (18a) with [<b,0] £ X we get 


a Ui • $ — / p 1 V $ = 


g • $. 

J n. 


The above holds for all <!> £ then, the equation (10a) follows in the H _1 (fli)-sense. Moreover, 

recalling that ui,g £ L 2 (fli), the strong constitutive Darcy equation (10a) also holds in the L 2 (Di)-sense. 
Now let ip £ C'g°(Di) and test equation (18a) with [0, gj] £ X to get — u 2 ' ^P = Jn 2 FP- Since this 
holds for all smooth functions the strong conservative equation (lOe) follows. 


Next we focus on the boundary and the interface conditions. Take vi £ Hdi v (^i), test (10a) with 
[vi, 0] £ X, integrate by parts and get 


n) dS. 


/ aui • vi + / p 2 (vi • n) dS + / Vpi • Vi - / p 1 (y 1 -V)dS=- g • Vi + / /s (vi • n) 

Jn i Jr J Oi Jon i Jth Jr 

We split the boundary term on in two pieces, the interfaces network T and the outer part d£l i — T = 
fl dft and replace u with n using the relationship given in Definition 11 equation (l8|). Additionally, 


recalling that the identity (10a) is satisfied, the expression above writes as 


[ P 2 (vi -n)dS - [ pi (vi -V)dS - f pi (vi ■ V) dS = [ / s (vi • n) dS. 

Jr Jr Jon^on Jr 

Since the above holds for all Vi £ Hdi v (^i) and the map Vj vi • n from Hdi v (^i) onto H~ 1 ^ 2 (dH i) is 
surjective, the normal stress balance condition across the interfaces (lOh) and the drained Dirichlet boundary 
condition (10c) follow in the sense of P 1 / 2 (r) and H 1 ^ 2 (dfliCidfl) respectively. Finally, taking q 2 £ H 1 (fl 2 ), 


testing (10a) with [ 0 , 52 ] € X and integrating by parts we get 


/ pp 2 q 2 dS— / (u 1 -n)q 2 dS+ V • u 2 q 2 - / (u 2 • v)q 2 dS = / F q 2 - l fft q 2 dS. 

Jr Jr Jn 2 Jon 2 Jn 2 Jr 


(25) 


/r Jr Jn 2 

Again, we split the boundary term on dfl 2 in two pieces, the interfaces network T and the outer part 


3D 2 —r = <9n 2 ri(9f2; then replace v with n using the relationship given in Definition 11 equation (|8j) . Since 
the identity (lOe I is satisfied, the expression above reduces to 

/ PP 2 Q 2 dS - / (U! • n) q 2 dS + / (u 2 • n)q 2 dS - (u 2 • n)q 2 dS = - / f A q 2 dS. 

Jr Jr Jr Jan 2 n n Jr 


(26) 


Observing that the above holds for all q 2 £ E( fl 2 ), it follows that the normal flux balance condition across 
the interfaces (lOg) and the null normal flux boundary condition (lOf) hold, in the sense of P -1 / 2 (T) and 
H~ 1 ^ 2 (dfl 2 D (9fl) respectively. This completes the proof. □ 


12 








































Finally in order to identify which forcing terms can be modeled using this formulation, this section closes 
characterizing the orthogonal projection onto the spaces V/f^) and V-^fD). 

Lemma 13. Let v £ L 2 (fl 2 ). Let £ £ Hq(Q) Q E(Sl 2 ) and ?y £ Eq(Q 2 ) be the unique solutions of the 
respective Dirichlet and Neumann problems 

— V ■ = — V • v in fl 2 , £ = 0 ondLl 2 . (27a) 


— V • V 77 = 0 infl 2 , V ?7 • n = (v — V£) • n ondfl 2 . (27b) 

Then, v — — V 77 is the projection of v onto V^f^) and + V?y is the projection of v onto 

PROOF. First recall that since C E(Q 2 ) then 

V _L (f2) = {v £ L 2 (D 2 ) : V • v = 0, in 1^2 and v • n = 0, on <9fI 2 }, 

i.e. V J -(fi 2 ) C Hdiv(f22). Next, observe that since £ £ iT^f^) then V£ £ L 2 (D 2 ), and because it is a 
solution to problem ( |27a[ ) it follows that V • (v — V£) = 0 i.e., v — € H div (fl 2 ). Moreover, for any 

q £ H 1 (£l 2 ) it holds that 


<(v- VO • n,q) 


R- i /2(an 2 ),ff-i/2( 9 n 2 ; 


= / (v-VO-Vg. 


In particular, if <7 = 1 then ((v — V^) • n, l) 1 


R- l /2(an 2 ),//-i/2(on 2 ) 


= 0 i.e., the data for the Neumann problem 


(27b) satisfy the compatibility condition and the problem has a unique solution 77 in Eq(Q 2 ). Now, it is 
clear that v -VZ-Vp £ H(n 2 ) and since v — (v — V£ — V?y) = V£ + V?y is orthogonal to H(fi2)> the 
result follows due to the characterization of orthogonal projections in Hilbert spaces. □ 


5. Conclusions and Final Discussion 


The present work yields several as summarized below. 


(i) The study of communication in complex networks, as stated in mm. can be perplexing due to the large 
number of nodes and links. One of the greatest achievements of this article is to provide a relatively 
simple upscaled description to an otherwise very complicated study. This result gives the grounds to 
use well-studied numerical methods such as Finite Elements on Complex Network Theory. 


(ii) The function a(-) presented in the strong form (10a), (10a), can be used as a scaling tool. If the volume 
of a region does not correctly reflect the impact of the universe it represents within the multiverse then, 
a(-) can be scaled in that region to counterbalance this deficiency and weigh each region of the map 
accordingly. 


(iii) The interface conditions (lOg) and (lOh), allow discontinuities along paths in the domain. This permits 
the modeling of sinks and/or sources on different manifolds along the region of interest. 


(iv) The method for PDE analysis on graphs presented in this article and the previous achievements in 
the literature are dual concepts. The preexisting results rely on complex definitions for the operators 
of the PDE and simple definitions for the domain. This work introduces simple definitions for the 
operators of the PDE and very restrictive definitions for the graph. As many results in mathematics 
have shown, complementing dual concepts and a wise interplay between them, are the key for a very 
strong theory. The result presented proves that this method for analyzing PDE on graphs is the 
gateway for significantly deeper understanding in the field. 
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(v) In the articles m and m, a similar mixed variational formulation is introduced in order to model the 
saturated fluid flow within a fractured porous medium. This result is an alternative version, allowing 
to treat fractures with substantially more general geometry. Additionally, the mixed variational for¬ 
mulation presented here is, to the authors’ best knowledge, unprecedented in the specialized literature 
of PDE analysis. 

For the limitations of the method we point out the following ones. 

(i) There are two particularly important cases of discretization for 3-D polygonal domains in a bipartite 
fashion, using tetrahedra as discussed in HE or using in cubes. For both cases, it is not difficult 
to prove that a subdivision of the graph /C 33 is contained in most of the graphs, defined by the 
“tetrahedral grids” or the “cubic grids”, associated to a given domain. Therefore, most of the time, 
the associated graph is nonplanar and bipartite. Although in these scenarios the domain for the PDE 
setting is already defined, the nonplanar structure of their natural associated graphs suggests that the 
analysis for bipartite, nonplanar graphs is an important issue still to be addressed. 

(ii) The variational formulation is dependent upon the domain chosen for the graph at hand. Thus, the 
construction of the domain that suits the application best will depend on the context. In this article 
we have suggested two possible choices, namely the tubular and downscaling maps. 

Finally, for further work we highlight the following aspects. 

(i) The method of analyzing highly clustered complex networks presented here is based upon the assump¬ 
tion that, approximating the behavior of a large number of particles with an open set, introduces an 
acceptable error when they are highly clustered. In order to understand the extent of this claim, fur¬ 
ther work about the error introduced by this approximation is of central importance. Notice that the 
quality of the approximation may very well depend on the nature of the quantities one is considering, 
which will impact in the PDE model. This aspect will be explored in future work. 

(ii) The Darcy Flow Model is a more general case of the Poisson equation whose unique differential operator 
is the Laplacian. Discrete versions of the Laplacian Operator for graphs are presented in [SUMS], a 
natural conjecture is that the discrete Laplacian operator will approximate the continuous version if 
the number of nodes is very large. This conjecture will be addressed in future work, because it may 
provide the foundation to justify the claim that, jumping from discrete to continuous is a reasonable 
estimate for the global behavior when the number of nodes is large. 

(iii) The formulation introduced divides the domain in two types of regions, namely fR and Q 2 which is 
the reason why the main setting of the problem had to be bipartite. Since every plane map is four 
coloreable it remains for future work to find a formulation addressing up to four types of regions. 
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